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We study the nucleon solution of the relativistic Faddeev equation as a function of 
density in the framework of a generalized Nambu— Jona-Lasinio model. We trun- 
cate the interacting two-body channels to the scalar diquark channel, the coupling 
constant of which is treated as a parameter. A 3-momentum cut-off is used to 
regularize the model. The Faddeev equation is solved numerically using the meth- 
ods developed by Tjon and others. At zero density the nucleon is bound only for 
unrealistically large values of the scalar coupling and this binding energy decreases 
quickly with increasing density. 



1 Introduction 

In contrast to high-temperature zero-density QCD, rather httle is known about 
high-density zero-temperature QCD. Due to technical difficulties (the fermionic 
determinant becoming complex at finite chemical potential), lattice calcula- 
tions are not able to provide unambiguous results. However, models of QCD 
seem to indicate a rich phase structure in high density quark matter. In par- 
ticular, n:iYp|h attention has recently been devoted to so-called colour supercon- 
ductivity Ol'BI: at high density, an arbitrarily weak attraction between quarks 
makes the quark Fermi sea unstable with respect to diquark formation and 
induces Cooper pairing of the quarks (diquark condensation). However, the 
groups who have studied colour superconductivity focused only on instabilities 
of the Fermi sea with respect to diquarks and have not considered possible 
3-quark clustering. To address this question would necessitate in principle a 
generalisation of the BCS treatment. As a first step we can look for insta- 
bilities of the quark Fermi sea with respect to 3-quark clustering by studying 
the evolution of the nucleon binding energy with density. A bound nucleon at 
finite density would be a signal of instability. In this study, we will use the 
Nambu- Jona-Lasinio (NJL) model and solve the relativistic Faddeev equation 
for the nucleon as a function of density. 
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2 The model 

The NJL model provides a simple implementation of dynamically broken chiral 
symmetry. It has been successful dp. the description of mesonic properties 
at low energy and several groups B'QB have used it to study baryons at zero 
density (for a review of the NJL model, see for example □) . Several versions 
of the NJL lagrangian are available. Whatever the version we choose, a Fierz 
transformation allows us to order the terms according to their symmetries in 
the qq channel; here we are only interested in the scalar and pseudoscalar 
terms: 

C^ = \g.[{i^i^f -{4>l5ri^f] ■ (1) 

To study the baryons, another Fierz transformation has to be performed in 
the qq channel. In this work we shall keep only the scalar diquark channel: 

Cs=gs{i>{i,C)T2p^i>'^){V{C-^l5)T2P^i^) , (2) 

where (3^ = 1/8/2 A'^ for A = 2,5,7 projects on the colour 3 channel and 
C = 17270 is the charge conjugation matrix. 

The ratio gs/ g-n depends on the version of the NJL lagrangian used. In 
the following, we will not choose a particular version of the model but rather 
leave the ratio gs/ g-n as a free parameter. We regularize the model with a 
3-momentum cut-off A. We have two parameters g.,, and A, which are fitted to 
the values of the pion decay constant = 93 MeV and the constituent quark 
mass M = 400 MeV. This gives us 5^ = 7.01 GeV-^ and A = 0.593 GeV. 

The effect of density is introduced by imposing the quark 3-momentum 
to be larger than the Fermi momentum kp. Solving the gap equation as a 
function oi kp gives us the usual dependence of the quark constituent mass on 
density. Chiral restoration occurs at kp/A — 0.58, which corresponds to about 
2.1 times the nuclear matter density 

3 Diquark at finite density 

As a first step to the resolution of the Faddeev equation, the scalar diquark 
mass has to be calculated as a function of density. The Bethe-Salpeter equation 
for the 2-body T-matrix is solved in the scalar qq channel using the ladder 
approximation (the explicit solution is given inSu) and the pole of the T-matrix 
gives the mass of the bound diquark. Note that the denominator of the scalar 
qq T-matrix is formally identical to that of the pionic T-matrix, except for the 
replacement of gs by g^^. That means that the pion and the scalar diquark are 
degenerate (and of zero mass) for a ratio gs/ g-n — 1- Fig- 1 gives the binding 
energy of the scalar diquark, Bdiq = 2y/kp + — Ediq as a function of the 
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Figure 1: Diquark binding energy (in GeV) as a function of kp/A 



dimensionless variable fci?/A for a value of the scalar coupling gs/g-K = 0.83. 
One can see that the scalar diquark is bound for all values of the density, the 
binding energy increasing with density to reach more than half a GeV. The 
sharp peak, which occurs at chiral symmetry restoration, is a consequence of 
the choice of a large scalar coupling, close to the value gs/g-ir — 1, for which 
the scalar diquark and the pion are degenerate. 



4 The relativistic Faddeev equation 

Because of the separability of the NJL interaction, the 3-body relativistic Fad- 
deev equation in the ladder approximation can be reduced to an effective 2- 
body Bethe-Salpeter equation describing the interaction between a quark and 
a diquark. Details concerning the derivation oLthis equation can be found in 
the previous studies performed at zero density era. Explicitly, it is written: 

*(^. I^^ij d\'R{2P/i + q')V{q, q'; P)^{P, q') , (3) 

where R{2P/3 + q') is the two-body T-matrix for the scalar diquark and 
V{q,q';P) involves the product of the propagators of the spectator and ex- 
changed quarks: 



yM-^n = P^^^^M^, (4) 



Here p'^ (i=l,2,3) arc the momenta of each valence quark, P — pi + p2 + ps is 
the total momentum of the nucleon, and g, q' are the Jacobi variables defined 
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by P3 = P/3 — q; p'l = P/3 — q' . Wc now look for the nuclcon solution of (H): 

\cr.q<i>2iqa,q) J ^ ' 

With this form for Eq. (|) becomes a set of two coupled integral equations. 
Following Huang and Tjona, we then perform a Wick rotation on the qo and 
q'g variables. Thia leads to two coupled complex integral equations which are 
solved iterativelyQ'Ll. The initial guess for each of the wave functions $i((7oj q) 
and $2 (90 1 9) consists of a Gaussian for the real part and a derivative of a 
Gaussian for the imaginary part. The number of iterations needed to reach 
convergence of the solutions is about four. 

The above discussion remains valid at finite density. Again, we incorporate 
the effects of density by restricting the 3-momentum of each valence quark to 
values larger than the Fermi momentum kp, i.e.: 

kF<\p^'\<^ (i = 1,2,3). (6) 

This condition translates into Eq. (|^) as a complicated cut-off on the 3-momen- 
tum integration variable. Apart from this restriction, the method of solving 
the Faddeev equation is the same as at zero density. 



5 Results and conclusions 

At zero density we found that the nucleon is bound only if the scalar coupling 
is strong enough, i.e. gs/d-n ^ 0-8, which means that the nucleon is not bound 
in either the "standard" NJL model [gs/g-n = 2/13) or the colour-current 
interaction version {gs/g-n = V^)- 

At finite density, we solve the Faddeev equation for the energy Enuc of the 
nucleon as a function of the Fermi momentum. Results are shown in Fig, [gl f or 
gs/gn — 0.83. The binding energy of the nucleon, Bnuc = Ediq + \/kp + AP — 
Enuc, is depicted, again as a function of the dimensionless variable kp/ K. The 
binding energy of the nucleon is relative to the quark-diquark threshold: as 
there is no confinement in the NJL model, nothing can prevent the existence of 
a free diquark. In contrast to the diquark (shown in Fig. ^ for the same value 
of gs/g-n), which is bound over the whole range of densities, the binding energy 
of the nucleon decreases quickly with density and the binding disappears well 
before nuclear matter density (which corresponds to kp/A ~ 0.45). 

These results imply that, in the region characteristic of colour supercon- 
ductivity (beyond chiral symmetry restoration), the quark Fermi sea is un- 
stable only with respect to formation of diquarks and not of 3-quark clusters. 
However, we have to emphasize that we included only the scalar part of the 
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Figure 2: Nucleon binding energy (in MeV) as a function of kp/K 



qq interaction; at zero density, several authors Su have shown that the axial- 
vector qq interaction gives an important contribution to the nucleon binding 
energy (of the order 100 MeV), while the axial- vector diquark is not bound for 
reasonable values of the axial-vector coupling. Incorporating the axial-vector 
qq interaction is necessary to obtain a more realistic picture of the evolution 
of the nucleon binding energy with density and could significantly modify the 
present results. 
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